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We study chaotic regions in the phase space of classical non-Abelian gauge theory, fo- 
cusing particularly on those which determine the low-energy interactions between BPS 
monopoles, and comment on the relevance of the obtained results for long-standing spec- 
ulations which relate classical Yang-Mills chaos to the disordered quantum vacuum and 
quark confinement. 
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, 1. Why study chaos in (semi-) classical gauge theories? 

Deterministic chaos in the time evolution of classical non-Abelian (and hence non- 



linear) gauge theories has been investigated for about thirty yeara^. This ongoing 
qq ■ endeavour has several motivations. On the more conceptual side, the chaotic be- 

havior of gauge theories reveals typical signatures of quantum chaos, visible e.g. in 
the distribution of nearest-neighbor level spacings of lattice Dirac spectra according 
to Gaussian matrix ensembles^, shows a continuous cascading of the dynamical de- 
grees of freedom (and their energy) towards the ultraviolet during time evolution^ 
and has potential implications for the continuum limit of lattice gauge theories^. 
One of the oldest and perhaps most important fundamental motivations for inves- 
tigating chaos in non-Abelian gauge theories, furt herm ore, was to shed light on its 
potential role in the QCD confinement mechanisirPE] Much of the work with more 
phenomenological focus, on the other hand, attempts to gain insight into otherwise 
hardly accessible non-equilibrium processes by exploiting the increasingly classical 
behavior of long-wavelength fields with growing temperature. Such processes are of 
great interest in the context of current experimental heavy-ion programs at RHIC 
and LHC which create and analyze matter under extreme conditions. The so far best 
studied examples are particle production from collective fields^ and relations be- 
tween the maximal Lyapunov exponents of t he c lassical dynamics and the damping 
or thermalization rates of hot gauge systemd^l 

Astrophysical and cosmological applications, motivated by the classical thermo- 
dynamics of hot and dense gauge theories as well, range from the description of 
stellar interiors to the role of chaos during semi-classical evolution phases of the 
early universe. Interesting examples are topological structure formation and its im- 
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pact on baryon number violating production processes in the Standard Model at 

[si 

temperatures around the electroweak phase transition . The chaoticity properties 
of particle motion in various curved spacetimes have also been investigated (see e.g. 
Ref.lU). 

A significant part of the results on gauge-theory chaos was obtained in Yang- 
Mills-Higgs (YMH) theories. These are well suited for the study of transitions from 
quantum to classical chaotic behavior since their weak-coupling and semiclassical 
limits are controllable at all length scales (in contrast to those of pure YM theories 
and QCD). The chaotic behavior of the YMH theory with Higgs fields in the fun- 
damental representation of the gauge group (which forms a part of the electroweak 
sector of the standard model) is e.g. of relevance for the mentioned baryon number 
violating processes. YMH theories with Higgs fields in the adjoint representation, on 
the other hand, appear in grand-unified theories and are of special interest because 
they sustain stable, finite-action monopole solutionis. 

The latter will be a main focus of this article. More specifically, we will study the 
regular and chaotic regimes in the low-energy dynamics of the two-monopole sys- 
tem, and the transitions between them. This dynamics is amenable to an enormous 
but controlled dimensional reduction of the relevant phase space - the geodesic 
approximation - which we will outline in the following section. In the subsequent 
Sec.|3]we quantify the chaoticity of the monopole- monopole interactions in various 
phase space regions, and in Sec.2]we will address potential implications for vacuum 
disorder and the quark confinement mechanism in quantum Yang-Mills theory and 
QCD. Section [5J finally, contains a summary and some conclusions. 



2. Geodesic two-monopole dynamics 

In the following two sections we are going to report on a recent investigation of 
regular and ch aotic low-energy interactions among monopoles in the simplest pos- 
sible settingpJ], i. e . between two electrically charged magnetic Bogomol'nyi-Prasad- 
Sommerfield (BPS) monopoles^l or dyontQ which solve the SU(2) YMH Bogo- 
mol'nyi equation^ 

Bt = = ± {S aC % + 9^° A) $ c (1) 

(where is the field strength tensor of the gauge field A® and $ a is the (adjoint) 
Higgs field). 

The two-BPS monopole system is prototypical for physically interesting sub- 
systems of the gauge dynamics whose spatially varying fields - here the solitonic 
monopoles with topologically induced magnetic charge - and time evolution can 
be studied without invoking either uncontrolled approximations (as e.g. the dras- 
tic "homogeneous approximation" used in the pioneering studies^) or requiring the 



a These solutions bear i ntere sting similarites to the BPS dyon constituents of caloron solutions 
with nontrivial holonomjEII which consist of a BPS monopole-antimonopole pair (for N c = 2). 
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elaborate lattice solution of the full, hyperbolic Yang-Mills-Higgs equation d^ l ^ l 
This is because the low-energy time evolution of the dyon pair is accurately de- 
scribed by the geodesic motion of a point on the manifold which the few collective 
degrees of freedom of the two-monopole solution sparP31, i.e. it is governed by ordi- 
nary differential equations. 

The basis for this geodesic approximation is that all two-monopole solutions form 
a family whose members are parametrized by continuous collective coordinates or 
"moduli" x a . (For the one-monopole solution, e.g., these are the three position 
coordinates of the center and an overall phase angle.) The corresponding moduli 
space M.2 is a manifold whose metric is induced by the metric on the space of all 
finite-energy field configurations, i.e. by the kinetic terms of the YMH Lagrangian, 
and known explicitly. Owing to energy conservation and the degeneracy of all static 
two-monopole solutions, the low-energy dynamics of two BPS dyons then describes 
geodesic motion of the associated point on Mi . After separating the center-of-mass 
motion and an overall phase (whose time dependence is associated with the total 
electric charge), the remaining four-dimensional internal part of the moduli 

space can be parametrized by three Euler angles (p and ip, which determine 
the orientation of the two-monopole system, and the distance variable g which 
measures (at large g) the separation between the two centers. The metric <?^ H ' on 

has been constructed explicitly by Atiyah and HitchiJlSl, and the resulting 
internal Lagrangian is that of a non-rigid body with distance-dependent "moments 
of inertia" around the body-fixed axea^, 

L AH (x, x) = y 5^ H) (x) x a x p = L AH (g, i>, g, 0, -0) , (2) 

where m is the reduced mass of the monopoles. Physically, the validity of the 
geodesic approximation implies that at small velocities (compared to the velocity 
of light) internal excitations (vibrations) and deexcitations (radiation) of the dyons 
can be neglected, i.e. the monopole pair adapts adiabatically to its interactions by 
deforming reversibly and scattering elastically. 

3. Regularity and chaos in monopole interactions 

The possibility for chaotic behavior depends on the number of integrals of the 
motion of the underlying dynamics, which must be less than the number of degrees 
of freedom. For the geodesic dynamics ([2]) of the two-monopole system with its four 
degrees of freedom, three independently conserved quantities are known explicitly. 
Those are the total angular momentum M 2 , the energy Eah and the generalized 
momentum p v conjugate to the coordinate <p which is cyclic, i.e. does not appear 
explicitly in the Lagrangian ([2]). For the two-monopole dynamics to be (Liouville) 
integrable would therefore require the existence of minimally one additional constant 
of the motion. Such a fourth conserved quantity indeed exists at least if the two 
monopoles remain infinitely separated, since then their individual electric charge 
is conserved. This situation changes when the two dyons begin to approach each 
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other and only their total charge remains conserved during Higgs-induced charge 
exchange. Then chaotic motion becomes possible and first numerical evidence for 
its existence in the two-dyon phase space was gathered in Refs. [TTJ 

In Ref. [11] we have then systematically analyzed regular and chaotic two- 
monopole motions on the basis of thirteen long-time phase space trajectories for 
which four-dimensional time series were generated by numerically integrating the 
equations of motion with high accuracy over typically 2 25 time steps. The initial 
data sets were chosen to cover a representative range of motion patterns and to 
explore the low-energy dyon interactions at different strengths. The resulting set of 
orbits includes sequences of trajectories whose decreasing minimal dyon separations 
interpolate between asymptotic dyon distances, where charge exchange becomes in- 
effective and the geodesic dynamics integrable, and relatively small minimal separa- 
tions for which the interactions are expected to become non-integrable. Hence these 
orbits allow to map out the order-chaos transition in the two-monopole system. 

The first part of our analysis consisted in constructing the Poincare sections of 
this orbit set from the Hamiltonian on the reduced four-dimensional phase space 
in which p v and A4 2 are conserved and act as fixed "external" parameters. Energy 
conservation then constrains all orbits to three-dimensional hypersurfaces. (Numer- 
ically, E, M 2 and p v were conserved up to deviations of order 10~ 12 .) Orbits with 
weak initial Coulomb attraction between the dyons cover a rather large range of g 
values with relatively moderate variations of the radial velocities which stay well 
inside the asymptotic region of approximately (or KAM-) integrable motion. Hence 
their Poincare sections (in the (g,p g ) plane at Mi = 0) consist of one-dimensional, 
continuous closed curves corresponding to quasiperiodic motions. When increas- 
ing the initial Coulomb attraction, the variations in dyon distance become smaller 
(tighter orbits) while their relative momenta vary more strongly. The increased at- 
traction also brings the dyons closer together. From a certain minimal inter-dyon 
distance £? min ~ 2n onward the section visibly spreads out into a broadly distributed 
scatter of points which eventually fills the (g,p g ) plane. This is a typical signature 
for the corresponding aperiodic orbit to have become chaotic. 

In order to investigate this chaoticity further, we have calculated high-resolution 
power spectra of the momentum conjugate to the dyon separation for selected or- 
bits. This allows for a more accurate distinction between quasiperiodic and chaotic 
motion patterns and yields quantitative information about the underlying scales. 
The resulting spectra indeed clearly separate quasiperiodic from irregular behav- 
ior. Two of the orbits, as expected those with the maximal initial Coulomb force 
between the dyons, were identified as chaotic. This substantially increases previous 
evidence that the relative low-energy motion of two BPS dyons admits, apart from 
the asymptotic g — > oo region, only three independent conserved quantities and 
turns out to be genuinely non-integrable. 

In addition, the power spectra characterize quasiperiodic dyon-pair orbits (i.e. 
those which remain close enough to the asymptotic region) quantitatively by es- 
tablishing the number of fundamental modes (two), determining their frequencies 
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and yielding the strength distribution over their various harmonics^. The restric- 
tion to the minimal number of quasiperiodic modes is rather widespread among 
nonlinear dynamical systems if they are sufficiently strongly coupled. The common 
expectation that nonlinear couplings between more than two fundamental modes 
increasingly turn quasiperiodicity into chaos may therefore apply to the two-dyon 
system as well and explain why we have encountered only two-mode-quasipcriodic 
and chaotic trajectories. 

In contrast to their rather complete specification of quasiperiodic motion pat- 
terns, power spectra do extract relatively little pertinent quantitative information 
from aperiodic orbits. Hence we have additionally calculated those quantities which 
perhaps most directly quantify the chaoticity of irregular motion patterns, i.e. the 
maximal Lyapunov exponents, for a suitable set of orbits. As expected, the Lya- 
punov exponents of orbits previously identified as quasiperiodic were found to van- 
ish. The two orbits with an irregular broadband power spectrum, on the other 
hand, turned out to have finite and positive maximal Lyapunov exponents whose 
values were approximately determined as L maXi i ~ 0.02 and L maXj2 ~ 0.008. These 
exponents provide our most unequivocal and quantitative evidence for the chaotic- 
ity of the dyon-dyon interactions. The orbit with the smaller Lyapunov exponent, 
furthermore, shows signs of intermittent behavioral 



4. Chaotic monopole systems and confinement 

In the following section we will summarize several ideas and speculations on the 
potential releva nce of the discovered two- monopole chaoticity for vacuum disorder 
and confinement! 1 8 | 1 9 ] m q uan tum YM and YMH theories. 

We start by recalling the long-standing conjecture that the vacuum of non- 
Abelian gauge theories, when undergoing a transition from weakly to strongly cou- 
pled fields, also undergoes an order-disorder transition, and that the strongly cou- 
pled QCD vacuum is populated by highly irregular color field configurations^. In the 
limit of a large number of colors, in particular, a vacuum made of random Yang- Mills 
fields is known to be a necessary and sufficient condition for quark confinement!^! 
From the outset, one of the motivations for investigating chaos in non-Abelian 
gauge the ories was therefore to shed light on its potential role in the confinement 
mechanisrrPJ^l More recently, a numerical investigation of the classical time evolu- 
tion of Yang-Mills field configurations generated by finite-temperature (quantum) 
lattice simulations^^ has provided evidence for the confining strong-coupling phase 
to be indeed substantially more chaotic than its weakly coupled counterpart: as a 
function of increasing coupling the maximal Lyapunov exponent undergoes a sharp 
transition to larger values around the confinement transition. 

Moreover, the instability of constant color-magnetic vacuum fields^l made it 
natural to assume that both gauge invariance and stability of the physical vacuum 
may be restored by disordering color-magnetic background fields. Under the gluonic 
infrared degrees of freedom^! envisioned to exhibit (and maybe cause) this disorder 
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are random domains as well as vacuum populations of quasi-randomly distributed, 

percolating center vortices or monopoles^. 

The pivotal role which non-Abelian monopoles and their chaotic interactions 

may play in the context of quark confinement in QCD is further presaged in TV = 2 

1041 

supersymmetric Yang-Mills theory in 3+1 dimensions^ where BPS monopoles in- 
deed realize a non-Abelian version of the classic 't Hooft-Mandelstam dual super- 
conductor confinement mechanism^! In this scenario the condensation of magnetic 
BPS monopole charges screens color magnetic charges and confines color electric 
charges by the dual Meifiner effect. Similar scenarios, in which the condensation of 
monopole-like objects plays a key role, are expected to unfold in more physical gauge 
theories as well. As a case in point, in the 2+1 dimensional Yang-Mill-Higgs model 't 
Hooft-Polyakov monopoles (a generalization of BPS monopoles which play the role 
of instantons in this case) generate "weak confinement" by forming a monopole an- 
timonopole plasma, as shown by Polyako\ESI. In 3+1 dimensional Yang-Mills theo- 
ries, furthermore, there is lattice evidence for the condensation of Abelian-projected 
monopoles to generate the bulk of the string tension ' . (According to an interesting 
recent suggestion, the "active" monopoles mi ght actually be BPS dyon constituents 
of caloron solutions with nontrivial holonomy^-SI) 

In light of the above arguments it is tempting to speculate that a disordered 
ensemble of monopoles (and anti-monopoles) in a semiclassical vacuum may be gen- 
erated by the classically chaotic low-energy interactions among individual monopole 
pairs which we have investigated above. Below we will suggest two approaches to- 
wards pursuing and testing such ideas in a more quantitative fashion. 



5. Summary and conclusions 

We have discussed chaotic regions in the classical phase space of non-Abelian gauge 
theories and studied, in particular, regular and chaotic motion patterns of two inter- 
acting BPS monopoles at low energies. Our analysis is based on a representative set 
of long-time phase-space trajectories in the geodesic approximation and intended to 
survey the order-chaos transition and to characterize the quasiperiodic and chaotic 
behavior quantitatively. 

The observed changes in the dimension of the trajectories' Poincare sections 
(from one to two) provide clear evidence for transitions from quasiperiodic to chaotic 
motion when the monopoles come close enough to each other. These results were 
confirmed and quantified by the analysis of high-resolution power spectra and Lya- 
punov exponents for selected orbits. The obtained values for the maximal Lyapunov 
exponents contain information on the relaxation time and thermalization (damp- 
ing) rate of a non-equilibrium dyon system at sufficiently high temperatures. Taken 
together, our results provide convincing evidence for and a quantitative description 
of both quasiperiodic and chaotic regions in the low-energy phase space of two BPS 
dyons. They also show that no more than the three explicitly known integrals of 
the motion are conserved by the geodesic forces between the dyons. 
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Since the motion of free dyons is integrable, the chaotic behavior analyzed above 
can be uniquely traced to the interactions between the dyons. This raises hopes that 
our quantitative understanding of the chaotic dyon-dyon interactions may also gen- 
erate new insights into the disorder of interacting monopole ensembles of the type 
which are expected to populate the vacuum of the strong interactions. For suffi- 
ciently dilute systems, expansions in the monopole density and more sophisticated 
many-body techniques may e.g. provide a quantitative treatment of chaotic multi- 
monopole systems. Another option would be the technically challenging extension of 
the geodesic approximation to approximate multi-monopole-antimonopole solutions 
with their more complex interactions. 
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